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The free energy of a field theory can be considered as a functional of the free correlation function. 
As such it obeys a nonlinear functional differential equation which can be turned into a recursion 
relation. This is solved order by order in the coupling constant to find all connected vacuum diagrams 
with their proper multiplicities. The procedure is applied to a multicomponent scalar field theory 
with a <^*-self-interaction and then to a theory of two scalar fields and A with an interaction (jy^A. 
All Feynman diagrams with external lines are obtained from functional derivatives of the connected 
vacuum diagrams with respect to the free correlation function. Finally, the recursive graphical 
construction is automatized by computer algebra with the help of a unique matrix notation for the 
Feynman diagrams. 



I. INTRODUCTION 



If one wants to draw all Feynman diagrams of higher orders by hand, it becomes increasingly difhcult to identify 
all topologically different connections between the vertices. To count the corresponding multiplicites is an even more 
tedious task. Fortunately, there exist now various convenient computer programs, for instance FeynArts or 
QGRAF |^,|^, for constructing and counting Feynman graphs in different field theories. These programs are based on 
a combinatorial enumeration of all possible ways of connecting vertices by lines according to Feynman's rules. 

The purpose of this paper is to develop an alternative systematic approach to construct all Feynman diagrams of 
a field theory. It relies on considering a Feynman diagram as a functional of its graphical elements, i.e. its lines and 
vertices. Functional derivatives with respect these elements are represented by graphical operations which remove 
lines or vertices of a Feynman diagram in all possible ways. With these operations, our approach proceeds in two 
steps. First the connected vacuum diagrams are constructed, together with their proper multiplicities, as solutions of 
a graphical recursion relation derived from a nonlinear functional differential equation. These relations have been set 
up a long time ago , but so far they have only been solved to all orders in the coupling strength in the trivial case 
of zero-dimensional quantum field theories. The present paper extends the previous work by developing an efficient 
graphical algorithm for solving this equation for two simple scalar field theories, a multicomponent scalar field theory 
with (/)^-self-interaction and a theory with two scalar fields (f> and A with the interaction (fPA. In a second step, all 
connected diagrams with external lines are obtained from functional derivatives of the connected vacuum diagrams 
with respect to the free correlation function. Finally, we automatize our construction method by computer algebra 
with the help of a unique matrix notation for Feynman diagrams. 

II. SCALAR <;/>*-THEORY 



Consider a self-interacting scalar field (f) with N components in d euclidean dimensions whose thermal fiuctuations 
are controlled by the energy functional 

E[4>]^\ I G^^(j)lC^2 + j: I Vi234M2<t'3^4 (2.1) 
^ J12 ^- J1234 

with some coupling constant g. In this short-hand notation, the spatial and tensorial arguments of the field 0, the 
bilocal kernel G^^, and the quartic interaction V are indicated by simple number indices, i.e., 

1 = {xi,Q;i}, y = X! Z'^''^;!, 01 = , G^2 ^ G~l,^^{xi,X2) , V1234 = ^ai, 02, 03,04 (2^1, 2:2, 2^3, 2^4) • (2.2) 

The kernel is a functional matrix G~^, while ^ is a functional tensor, both being symmetric in their indices. The 
energy functional (2.1) describes generically c?-dimensional euclidean 0*-theories. These are models for a family of 
universality classes of continuous phase transitions, such as the 0(iV)-symmetric 0^-theory which serves to derive 
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the critical phenomena in dilute polymer solutions (A'' = 0), Ising- and Heisenberg-like magnets {N = 1,3), and 



superfluids {N = 2). In all these cases, the energy functional (2J_) is specified by 

Gala2i^U^2) ^ Sa,,a2 {~d'^^ + m'^) S{xi - X2) , (2.3) 
,a2 ,CK3 ,<^4 ("^l 1 "^2 , -^3 ; -^4) — ^ {^cti ,a2 ,Q4 ~t~ (5q^ (5q2 ,Q4 ~t~ ^ai ,a4 <^Q2 ,a3 } ^{-^1 -^2 )<^(^1 -^3 )^('^1 -^4 ) , (^■^) 

where the mass is proportional to the temperature distance from the critical point. In the following we shall leave 
and V co mpletely general, except for the symmetry with respect to their indices, and insert the physical values 



(2.3) and (2.4) at the end. By using natural units in which the Boltzmann constant kg times the temperature T 
equals unity, the partition function is determined as a functional integral over the Boltzmann weight e^^I"^' 

P^e-^I-^l (2.5) 

and may be evaluated perturbatively as a power series in the coupling constant g. From this we obtain the negative 
free energy = In Z as an expansion 

^-E^(lf) (2-6) 

The coefficients W'^p'^ may be displayed as connected vacuum diagrams constructed from lines and vertices. Each line 
represents a free correlation function 

1 2 = Gi2 , (2.7) 



which is the functional inverse of the kernel G ^ in the energy functional (2.1), defined by 

^012623! =<5i3. (2.8) 
The vertices represent an integral over the interaction 



1^234 . (2.9) 



1234 



To construct all connected vacuum diagrams contributing to PF^p) to each order p in perturbation theory, one connects 
p vertices with Ap legs in all possible ways according to Feynman's rules which follow from Wick's expansion of 
correlation functions into a sum of all pair contractions. This yields an increasing number of Fcynman diagrams, each 
with a certain multiplicity which follows from combinatorics. In total there are Al^pl ways of ordering the 4p legs of the p 
vertices. This number is reduced by permutations of the legs and the vertices which leave a vacuum diagram invariant. 
Denoting the number of self-, double, triple and fourfold connections with S,D,T,F, there are 2!'^, 2!^, 3!"^, 4!^ leg 
permutations. An additional reduction arises from the number TV of identical vertex permutations where the vertices 
remain attached to the lines emerging from them in the same way as before. The resulting multiplicity of a connected 
vacuum diagram in the (/)^-theory is therefore given by the formula 



The superscript E records that the number of external legs of the connected vacuum diagrams is zero. The diagram- 



matic representation of the coefficients W'^p^ in the expansion (2.6) of the negative free energy W is displayed in Table 
lup to five loops pOMTl. 



For higher orders, the factorially increasing number of diagrams makes it more and more difficult to construct all 
topologically different diagrams and to count their multiplicities. In particular, it becomes quite hard to identify by 
inspection the number N of identical vertex permutations. This identification problem is solved by introducing a 
uniqued matrix notation for the graphs, to be explained in detail in Section IV. 

In the following, we shall generate iteratively all connected vacuum diagrams. We start in Subsection II. A by iden- 
tifying graphical operations associated with functional derivatives with respect to the kernel G~^, or the propagator 
G. In Subsection II. B we show that these operations can be applied to the one-loop contribution of the free partition 



function to generate all perturbative contributions to the partition function (2.5). In Subsection II. C we derive a 
nonlinear functional differential equation for the negative free energy W, whose graphical solution in Subsection II. D. 
yields all connected vacuum diagrams order by order in the coupling strength. 
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A. Basic Graphical Operations 



Each Feynman diagram is composed of integrals over products of free correlation functions G and may thus be 
considered as a functional of the kernel . The connected vacuum diagrams satisfy a certain functional differential 
equation, from which they will be constructed recursively. This will be done by a graphical procedure, for which we 
set up the necessary graphical rules in this subsection. First we observe that functional derivatives with respect to the 
kernel or to the free propagator G correspond to the graphical prescriptions of cutting or of removing a single 
line of a diagram in all possible ways, respectively. 



1. Cutting Lines 



Since </> is a real scalar field, the kernel G "'^ is a symmetric functional matrix. This property has to be taken into 
account when performing functional derivatives with respect to the kernel G~^, whose basic rule is 



5G^l 



1 



{<5l3<542 + 5l4(532} 



(2.11) 



^From the identity (p.8|) and the functional chain rule, wc find the effect of this derivative on the free propagator 

(2.12) 



6Gi2 

2 _ _^ — G13G42 + G14G32 



5G 



34 



This has the graphical representation 
„ 5 



5G^l 



2 = 



(2.13) 



Thus differentiating a propagator with respect to the kernel G^^ amounts to cutting the associated line into two 
pieces. The differentiation rule ( p.ll ) ensures that the spatial indices of the kernel are symmetrically attached to the 
newly created line ends in the two possible ways. When differentiating a general Feynman integral with respect to 
G~^, the product rule of functional differentiation leads to a sum of diagrams in which each line is cut once. 

With this graphical operation, the product of two fields can be rewritten as a derivative of the energy functional 
with respect to the kernel 



»192 



= 2 



6E[(I)\ 

Jg^' 



(2.14) 



as follows directly from (2.1) and (2.11). Applying the substitution rule (2.14) to the functional integral for the fully 
interacting two-point function 



we obtain the fundamental identity 



Z 



G12 



5W 



(2.15) 



(2.16) 



Thus, by cutting a line of the connected vacuum diagrams in all possible ways, we obtain all diagrams of the fully 
interacting two-point function. Analytically this has a Taylor series expansion in powers of the coupling constant g 
similar to (2.6) 



G-12 — 



°° 1 / \ P 



p=o ■ 



(2.17) 



with coefficients 
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(2.18) 



The cutting prescription ( ^.16 ) converts the vacuum diagrams of pth order in the coefficients W^f-* in Table | to 
the corresponding ones in the coefficients G^j'' of the two-point function. The results are shown in Table || up to 
four loops. The numbering of diagrams used in Table || reveals from which connected vacuum diagrams they are 
obtained by cutting a line. For instance, the diagrams #15.1-#15.5 and their multiplicities in Table ^ follow from 
the connected vacuum diag ram # 15 in Table |. We observe that the multiplicity of a diagram of a two-point function 
obeys a formula similar to (2.1C): 



M. 



E=2 



4!Pp!2! 



(2.19) 



In the numerator, the 4:\p p\ permutations of the 4p legs of the p vertices are multiplied by a factor 2! for the permu- 
tations of the E = 2 end points of the two-point function. The number TV in the denominator counts the identical 
permutations of both the p vertices and the two end points. 



Performing a differentiation of the two-point function (2.15) with respect to the kernel G ^ yields 

SGi2 



- 2 



Gi234 — G12G3 



where G1234 denotes the fully interacting four-point function 



G 



1234 



1 
Z 



P0(/.l(/.20304e-^[*l 



(2.20) 



(2.21) 



The term G12G34 in (2.2C) subtracts a certain set of disconnected diagrams from G1234. By subtracting all discon- 
nected diagrams from G1234, we obtain the connected four-point function 



^1234 = G1234 — G12G34 — G13G24 ~ G14G23 



(2.22) 



in the form 



G 



1234 



6G1 

scT 



G13G24 — G14G23 . 



(2.23) 



34 



The first term contains all diagrams obtained by cutting a line in the diagrams of the two-point-function G12. The 
second and third terms remove from these the disconnected diagrams. In this way we obtain the perturbative expansion 



G 



1234 



p=i ^ 



4! 



ic,(p) 
''1234 



(2.24) 



with coefficients 



C,(P) 



234 



= -2 



SG^ 

SG ^ 



34 



p 

E 

9=0 



"'13 "'24 



/-<(p-9)/^(9) 
"14 "23 



(2.25) 



They are listed diagrammatically in Table III up to three loops. As before in Table II, the multiple numbering in 
Table [II indicates the origin of each diagram of the connected four-point function. For instance, the diagram #11.2.2, 
#11.4.3, #14.1.2^ #14.3.3 in Table |l| stems together with its multiplicity from the diagrams #11.2, #11.4, #14.1, 
#14.3 in Table p. 



The multiplicity of each diagram of a connected four-point function obeys a formula similar to (2.19): 

4!Pp!4! 



M. 



E=4 



(2.26) 



This multiplicity decomposes into equal parts if the spatial indices 1, 2, 3, 4 are assigned to the E — A end points of 
the connected four-point function, for instance: 
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62208 XXX - 20736 J XXX ^ +^0736 ^ XXX ' +^0736 ^XXX3 " (2.27) 



Generalizing the multiplicities (2.10), (2.19), and (2.26) for connected vacuum diagrams, two- and four-point fimctions 



to an arbitrary connected correlation function with an even number E of end points, we see that 

where N counts the number of permutations of vertices and external lines which leave the diagram unchanged. 

2. Removing Lines 

We now study th e gra phical effect of functional derivatives with respect to the free propagator G, where the basic 

= i {'5l3'542 + 5l4<532} . (2.29) 



differentiation rule ( 2.11 ) becomes 

6Gu ^ 1 
SGsi " 2 



We represent this graphically by extending the elements of Feynman diagrams by an open dot with two labeled line 
ends representing the delta function: 

1 — 2 6i2. (2.30) 



Thus we can write the differentiation ( 2.29 ) graphically as follows: 

1 2 ^XSl^S 4^2 + 1^4 3^2>. (2.31) 

^3 4 2 1 J 

Differentiating a line with respect to the free correlation function removes the line, leaving in a symmetrized way the 
spatial indices of the free correlation function on the vertices to which the line was connected. 

The effect of this derivative is illustrated by studying the diagrammatic effect of the operator 

[ Gi2-^. (2.32) 

Jl2 <JUi2 

Applying L to a connected vacuum diagram in W^p\ the functional derivative S/SG12 removes successively each of 
its 2p lines. Subsequently, the removed lines are again reinserted, so that the connected vacuum diagrams W^^^ are 
eigenf unctions of L, whose eigenvalues 2p count the lines of the diagrams: 

LW^P^ =2pW'^PK (2.33) 

As an example, take the explicit first-order expression for the vacuum diagrams, i.e. 

W'-'^=3[ 1/1234 G12G34, (2.34) 

J 1234 



and apply the basic rule ( p. 291 ), leading to the desired eigenvalue 2. 



B. Perturbation Theory 

Field the oret ic perturbation expressions are usually derived by introducing an external current J into the energy 



functional (2^) which is linearly coupled to the field (j) |13Hl5|. Thus the partition function ( |2.5| ) becomes in the 
presence of J the generating functional Z[J] which allows us to find all free n-point functions from functional derivatives 
with respect to this external current J. In the normal phase of a (/)^-theory, the expectation value of the field cj) is 
zero and only correlation functions of an even number of fields are nonzero. To calculate all of these, it is possible 
to substitute two functional derivatives with respect to the current J by one functional derivative with respect to the 
kernel G~^. This reduces the number of functional derivatives in each order of perturbation theory by one half and 
has the additional advantage that the introduction of the current J becomes superfluous. 
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1. Current Approach 



Recall briefly the standard perturbative treatment, in which the energy functional (2J) is artificially extended by 
a source term 



The functional integral for the generating functional 

Z[J] = J X>(/)e-^['^"^l 
is first explicitly calculated for a vanishing coupling constant yielding 

ZW[J] =exp|-iTrlnG-i + i^^ G12 JiJsj , 
where the trace of the logarithm of the kernel is defined by the series |16, p. 16] 



TrlnG-i = ^ / {G^,' - 5u} ■ ■ ■ {G~^ - 

„_i Jl...n 



(2.35) 



(2.36) 



(2.37) 



(2.38) 



If the coupling constant g does not vanish, one expands the generating functional Z[J] in powers of the quartic 
interaction V, and reexpresses the resulting powers of the field within the fun ctional integral (|]3£) as functional 
derivatives with respect to the current J. The original partition function ( ^.5[ ) can thus be obtained from the free 
generating functional ( 2.37] ) by the formula 



234 ^^^^ 6Ji5J25Ji6Ji 



(2.39) 



,7=0 



Expanding the exponential in a power series, we arrive at the perturbation expansion 

5^ 



1 



-9 



4! 
2 I 4! 



1234 ^^^^'' SJ15J2SJ35J4 

^ ^ J8 

12345678 ^^^'^ ^""^^ 5Ji5J25J'i5J45J;^5jQ5Jj5J9, 



(2.40) 



.7=0 



in which the pth order contribution for the partition function requires the evaluation of Ap functional derivatives with 
respect to the current J. 



2. Kernel Approach 

The derivation of the perturbation expansion simplifies, if we use functional derivatives with respect to the kernel 
in the energy functional (2T) rather than with respect to the current J. This allows us to substitute the previous 
expression (2.39) for the partition function by 



exp^-|/ F1.34 ,^-1,^-1 , (2.41) 

J1234 0G^2 0^34 



where the zeroth order of the negative free energy has the diagrammatic representation 

^'(O) = -^TrlnG-i = i Q . (2.42) 
Expanding again the exponential in a power series, we obtain 

^={' + ^/234 '^""W^^K^) /2345ar8 ^G-^G3-^G-^G- + • • ■ • ^'-"'^ 
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Thus we need only half as many functional derivatives than in ( ^.40 ). Taking into account (2.11). (2.12), and (2.38) 
we obtain 



1 



- = - Gl2 , 



1 



- — T {613024 + G14G23} 



such that the partition function Z becomes 



(2.44) 



^1234^5678 



12345678 



X 



9 ^12034056 078 + 24 O15O26O37O48 + 72 O12O35O46O78 



This has the diagrammatic representation 

2 



Z=<1 



4! 



4! 



oQOCO 



24 



72 OX) 



(2.45) 



(2.46) 



AH d iagrams in this expansion follow directly by successively cutting lines of the basic one-loop vacuum diagram 
( ^.42 ) according to ( ^.43 ). By going to the logarithm of the partition function Z , we find a diagrammatic expansion 
for the negative free energy W 



W 



4! 



4! 



24 



72 



(2.47) 



which turns out to contain precisely all connected diagrams in ( |2.46| ) with the same multiplicities. In the next section 
we show that this diagrammatic expansion for the negative free energy can be derived more efficiently by solving a 
differential equation. 



C. Functional Differential Equation for W = \nZ 



Regarding the partition function Z as a functional of the kernel O ^, we derive a functional differential equation 
for Z. We start with the trivial identity 



0, 



(2.48) 



which follows via functional integration by parts from the vanishing of the exponential at infinite fields. Taking into 
account the explicit form of the energy functional (|]^), we perform the functional derivative with respect to the field 
and obtain 



I i?^|<5i2-^or3\ 



6 



^1345 </'2</'3'/'4</'5U-^["*l -0. 



345 



(2.49) 



Applying the substitution rule (2.14), this equation can be expressed in terms of the partition function (2.5) and its 
derivatives with respect to the kernel O^^: 



S,2Z + 2 I Gll^^\g 
dG23 



Vl345 



S'Z 



345 



(2.50) 



Note that this linear functional differential equation for the partition function Z is, indeed, solved by (2.41) due to 
the commutation relation 



78 



^5678 TT^n exp 



(50. 



78 



Vi 



234 



1234 



SG^iSG^l 



(2.51) 
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which follows from the canonical one 

5 



Gr,} — G 



1 



{S13S24: + S14S23} 



Going over from Z toW = \nZ, the linear functional differential equation ( 2.50 ) turns into a nonlinear one: 

5W 5W 



^''^^ I ^''SG^~3 ' 7345 '''' I ^^^^ ^SG^SG^ 



'23 '^'-'45 



'23 '^'-'45 



(2.52) 



(2.53) 



If the coupling consta nt g vanishes, this is immediat ely s olved by ( 2.42| ). For a non- vanishing coupling constant g, 
the right-hand side in ( 2.53| ) produces corrections to (2.42) which we shall denote with Thus the negative free 

energy W decomposes according to 



(2.54) 



Inserting this into ( 2.53 ) and taking into account ( ^.44 ), we obtain the following functional differential equation for 
the interaction negative free energy 



G-, 



1 SW^ 



int 



SG12 4 J1234 

, 9 



1^234 G12G 



34 



/ 

J 1234 



1^234 G12 



~SG^ 



J2yy(int) JVK^'"*^ (5T/F('"*^ 

3 J^^ ■'1234 ^ x/^-i xr^-i ~^ x/^-i 



Vr. 



234 



SG12 SG^^ SG^2 f^G^ 



(2.55) 



With the help of the functional chain rule, the first and second derivatives with respect to the kernel G ^ are rewritten 
as 



r^-l tT13(-'24 77;— 

0(^12 J34 0tj34 



(2.56) 



and 



SG 12 ^G^l J 5678 



Gi5G26G37G'48- 



5^ 



SG^eSGjs 

7: I {G13G25G4Q + G14G25G3Q + G23G15G46 + G24G15G36} 



5G, 



56 



(2.57) 



respectively, so that the functional differential equation ( ^.55 ) for lI/^('"') takes the form (compare Eq. (51) in Ref. 0|) 

12 **Jl234 



G12 



12 



5G 



"T / ^1234 G12G34 — g / V1234 Gi2G35G46^-=; 

4 J1234 J123456 "^56 



„ , V1234 G15G26G37G48 < -r-; j-^ 1 j-^ } ■ (2.58) 

•J J 12345678 |^ 0Ct560Ct78 OLt^q 0Cr78 
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D. Recursion Relation and Graphical Solution 



We now convert the functional differential equation ( [2.58 ) into a recursion relation by expanding VF^'"*-' into a 
power series in G: 



p—i ^ 



(2.59) 



Using the property ( 2.33| ) that the coefficient satisfies the eigenvalue problem of the line numbering operator 

(2.32), we obtain the recursion relation 
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^(P+1)^12/ T/l234 Gl2G35G 

J123456 
p-1 



46^7=; / ^1234 G15G26G37G48T7; 77; — 

0(^56 J12345678 0Lr56''Lf78 



9=1 



P 

1 J 12345678 



14.234 G15G26G37G48 ■ 



<5G56 5G' 



78 



(2.60) 



and the initial condition (2.34). With the help of the graphical rules of Subsection II A, the recursion relation ( 2.60| ) 
can be written diagrammatically as follows 



61 2 (53 4 



+ 12 



5i — 2 



4 E' ^ 



q=l 



5i — 2 



1 X/ 3 ,5W^(«) 

2--^4 5Z 4 



This is iterated starting from 



p > 1, 



(2.61) 



= 3 



(2.62) 



The right-hand side of ( 2.61 ) contains three different graphical operations. The first two are linear and involve one 
or two line amputations of the previous perturbative order. The third operation is nonlinear and mixes two different 
one-line amputations of lower orders. 



An alternative way of formulating the above recursion relation may be based on the graphical rules 

1 




P 



SG12 




P 



(5Gi2(5G34 




P 



(2.63) 



With these, the recursion relation (2.61) reads 




P+l] = 4 
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p-i 

q=l 




(p-q 



p>i. 



(2.64) 



To demonstrate the working of (2.61), we calculate the connected vacuum diagrams up to five loops. Applying the 
linear operations to (2.6C), we obtain immediately 



Si 2 Si 2 S3 4 

Inserted into ( ^.61 ), these lead to the three-loop vacuum diagrams 



(2.65) 



W 



(2) _ 



24 



72 OX) 



(2.66) 



Proceeding to the next order, we have to perform one- and two- line amputations on the vacuum graphs in ( 2.66 ) 
leading to 



Si — 



= 96 



2 +144 



144 



(2.67) 
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and subsequently to 



Si — 2 53- 



288 



144 



288 2- 



144 3- 



144 1 



4 + 144 



(2.68) 



Inserting (2.67) and (2.68) into (2.61) and taking into account (2.65), we find the connected vacuu m di agrams of 
order p — 3 with their muhiphcities as shown in Table ^. We observe that the nonlinear operation in (2.61) docs not 
lead to topologically new diagrams. It only corrects the multiplicities of the diagrams generated from the first two 
operations. This is true also in higher orders. The connected vacuum diagrams of the subsequent order p = 4 and 
their multiplicities are listed in Table |. 

As a cross-check we can also determine the total multiplicities M^p'' of all connected vacuum diagrams contributing 
to W''P\ To this end we recall that each of the M^^^ diagrams in consists of 2p lines. The amputation of one or 
two lines therefore leads to 2pM(p) and 2p{2p- 1)M'-p^ diagram s wit h 2p—l and 2p—2 lines, respectively. Considering 
only the total multiplicities, the graphical recursion relations ( ^.6l| ) reduce to the form derived before in Ref. 



M^P+^) = 16p(p + l)Af(p) + 16^ 



These are solved starting with the initial value 



l)!(g-l)! 



p>l 



leading to the total multiplicities 

m(2) = 96 , m(3) = 9504 , M^"*) = 1880064 , 
which agree with the results listed in Table |. In addition we note that the next orders would contain 

M(^) = 616108032 , Af = 301093355520 , Af = 205062331760640 
connected vacuum diagrams. 



(2.69) 



(2.70) 



(2.71) 



(2.72) 



III. SCALAR 0M-THEORY 

For the sake of generality, let us also study the situation where the quartic interaction of the ^''-theory is generated 
by a scalar field A from a cubic (/>^A-interaction. The associated energy functional 



decomposes into the free part 



and the interaction 



i/ 1 2 1 2 

^ J 123 



Indeed, as the field A appears only quadratically in (3.1), the functional integral for the partition function 



(3.1) 



(3.2) 



(3.3) 



(3.4) 
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can be exactly evaluated with respect to the field A, yielding 

Z = y'pe^^'"'[^l (3.5) 

with the effective energy functional 

S^'^'^M = -^Trlni/-! + i / G^a' (/.i^2 - f / V125V346H56 M2h<l>4 ■ (3.6) 

^ ^ J12 O J123456 

Apart fro m a trivial shift due to the negative free energy of the field A, the effective energy functional (|]^) coincides 
with one ( ^.l[ ) of a (/)'*-theory with the quartic interaction 

^1234 = -3 / Vi25V3i6H5e. (3.7) 
J 56 



If we supplement the previous Feynman rules (2.7), (2.9) by the free correlation function of the field A 

1 2 = Hi2 (3.8) 

and the cubic interaction 



Vi23 , (3.9) 



123 



the intimate relation (3.7) between the (^''-theory and the (p'^A- theory can be graphically illustrated by 




(3.10) 

This corresponds to a photon exchange in the so-called s-, t- and w-channels of Mandelstam's theory of the scattering 
matrix. Their infinite repetitions yield the relevant forces in the Hartree, Fock and Bogoliubov approximations of 
many-body physics. In the following we analyze the (fi^ A-theory along similar lines as before the (/i^-theory. 

A. Perturbation Theory 

Expanding the exponential in the partition function ( ^.4| ) in powers of the coupling constant g, the resulting 
perturbation series reads 

^ = E TTli (!) / ^-^^^ / V,23V,56M2M5A3Ae e-^'"'[^-^l . (3.11) 



^0 (^py- ^4 



123456 



Substituting th e pro duct of two fields ip ot A by a functional derivative with respect to the kernels G ^ oi H ^, we 
conclude from ( 3.11 ) 

where the zeroth order of the negative free energy reads 

= -iTrlnG-i - IttIuH-' ^ i Q + ^ O ' (^'l^) 



Inserting (3.13) in (3.12), the first-order contribution to the negative free energy yields 



VF(1)=2 / Vi23V456 i?36Gi4G25+ / 1^123 ■^^456 i?36Gi2G45 , (3.14) 

J123456 J123456 

iagrams 

VF(i)=2 + r} — (~) . (3.15) 



'123456 J 123456 

which corresponds to the Feynman diagrams 
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B. Functional Differential Equation for W — \nZ 



The derivation of a functional differential equation for the negative free energy W requires the combination of two 
independent steps. Consider first the identity 



V 



which immediately yields with the energy functional (3.1) 



S{{A,)Z} 



'134 - 0, 

/3 OLT23 J34 0Lf23 



(3.16) 



(3.17) 



where (A) denotes the expectation value of the field A. In order to close the functional differential equation, we 
consider the second identity 



(3.18) 



which leads to 



{A,)Z = ^ [ V23iH,^4^. (3.19) 

J234 "^23 

Inserting (III) in (I3I7I ), we result in the desired functional differential equation for the negative free energy W = In Z: 

(3.20) 



612 + 2 I Gi3 -—J- - -2g / 1^134^467^^47 <^ .^-1.^-1 + 77^77^ 



23 



A subsequent separation ( ^.54 ) of the zeroth order ( [3.13 ) leads to a functional differential equation for the interaction 

na.rt nf thp frpp priprcrv • 



part of the free energy W 



/ ^12 ^ ~T / Vi23V456-ff36 {^12645 + 2G'i4G25} + 3 / V123V45QG12H; 

Jl2 0G12 * J123456 ^123456 



^123^56-^36 



123456 



123456 



(5G'45^ 



(3.21) 



Taking into account the functional chain rules ( ^.56 ), (2.57), the functional derivatives with respect to G ^ in (3.21) 
can be rewritten in terms of G: 



12^'" '5Gi2 4 



/ ^123^56-^36 {G12G45 + 2G14G25} + (7 / Vi23V456-ff36 {G12G47G58 

J123456 J123456 



-2G14G27G58} 



SG 



78 



Vi23V456i?36G'i7G28G39G, 



41 



1234567891 



J2^(int) JP7(int) SW'-™^'> 



(3.22) 



C. Recursion Relation And Graphical Solution 



The functional differential equation ( p.22| ) is now solved by the power series 

00 

^ (2p)! U; 

p— 1 ^ ^ ' 



(3.23) 



Using the property (2.33) that the coefficients H/^'') satisfy the eigenvalue condition of the operator ( 2.32| ), we obtain 
both the recursion relation 
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= 4(2p + 1) ( / V123V456 i?36 (G12G47G58 + 2G14G27G58) 

IJ12345678 



SGjs 



1234567891 



V V H r r r r ( J^W^ (2p\ SW^^-^) 5Wi^)\ 

''123 •'456 -n36LTi7LT28(-f39(-f41 | TFi TFi r / , 



SGjsSGgi ^ \2g / SGys SGqi J 



(3.24) 



and the initial value (3.14). Using the Feynman rules (2.7), (3.8) and (3.9), the recursion relation (3.24) reads 
graphically 



=4(2p+l)< 



SW^P^ 1 
Si — 2 2 



61 2 „ 



+ 



5l 2(53 4 3- 

4 ■ 



j2 f'^p\ SW'^P-i^ 2 

9=1 



2q 



Si 2 1 



4 (5VF(«) 

3 ^3 4 



p>l, (3.25) 



which is iterated starting from ( 3.15| ). In anal ogy t o ( 2.64 ), this recursion relation can be cast in a closed diagrammatic 
way by using the alternative graphical rules (2.63): 



p + 1 = 4(2p+ 1 






P > 1 



(3.26) 



We illustrate the procedure of solving the recursion relation (3.25) by construction the three-loop vacuum diagrams. 
Applying one or two functional derivatives to (3.15), we have 



_6W^ 

Si — 2 



Si 2 S3- 



(3.27) 



This is inserted into ( 3.25| ) to yield the three- loop diagrams shown in Table IV with their multiplicities. The table 
also contains the subsequent 4-loop results which we shall not derive here in detail. Obs erve that the multiplicity of 
a connected vacuum diagram D in the 0^yl-theory is given by a formula similar to (2.10) in the (/)^-theory: 



M 



(2p)!4P 



(3.28) 



Here 5* and D denote the number of self- and double connections, whereas N represents the number of identical vertex 
permutations. 



The connected vacuum diagrams of the ^^A-theory in Table IV can, of course, be converted to corresponding 
ones of the 0^-theory in Table ||, by shrinking wiggly lines to a point and dividing the resulting multiplicity by 3 in 
accordance with ( 3.10 ). This relation between connected vacuum diagrams in 0'*- and cjPA-thcory is emphasized by 
the numbering used in Table |rv|. For instance, the shrinking converts the five diagrams #4.1-#4.5 in Table IV to the 
diagram #4 in Table ||. Taking into account the different combinatorial factors in the expansion (2_^) and ( 3.23 ) as 
well as the factor 3 in the shrinkage (3.10), the multiplicity M^4^^ of a (^''-diagram results from the corresponding 



one M^2~A of the (/i^A-partner diagrams via the rule 



M. 



,fE=0 



1 



(2p-l)!! 



(3.29) 
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IV. COMPUTER GENERATION OF DIAGRAMS 



Continuing the solution of the graphical recursion relations (2.61) and (3.25) to higher loops becomes an arduous 
task. We therefore automatize the procedure by computer algebra. Here we restrict ourselves to the (/)^-theory because 
of its relevance for critical phenomena. 



A. Matrix Representation of Diagrams 

To implement the procedure on a computer we must represent Feynman diagrams in the ^^-theory by algebraic 
symbols. For this we use matrices as defined in Refs. |^,^. Let p be the number of vertices of a given diagram and 
label them by indices from 1 to p. Set up a matrix M whose elements Mij (0 < i,j < p) specify the number of lines 
joining the vertices i and j. The diagonal elements Ma {i > 0) count the number of self-connections of the ith vertex. 
External lines of a diagram are labeled as if they were connected to a single additional dummy vertex with number 
0. As matrix element Mqq is set to zero as a convention. The off-diagonal elements lie in the interval < My < 4, 
while the diagonal elements for i > are restricted by < Ma < 2. We observe that the sum of the matrix elements 
Mij in each row or column equals 4, where the diagonal elements count twice, 



j=0 



M,,, 



i=0 



A%-4. 



(4.1) 



The matrix M is symmetric and is thus specified by { p+ l) {p + 2)/2— 1 elements. Each matrix characterizes a unique 
diagram and determines its multiplicity via formula ( 2.28| ). From the matrix M we read offdirectly the number of 
self-, double, triple, fourfold connections S,D,T,F and the number of external legs E = J^^^i Moi- It also permits 
us to calculate the number N of identical vertex permutations. For this we observe that the matrix M is not unique, 
since so far the vertex numbering is arbitrary. In fact, N is the number of permutations of vertices and external lines 
which leave the matrix M unchanged [compare to the statement after ( 2.28| )]. If um denotes the number of different 
matrices representing the same diagram, the number N is given by 



r) 71 ,r J- J- 



(4.2) 



The matrix elements Moi count the number of external legs connected to the ith vertex. Inserting Eq. ( |4.2| ) into 
the formula (2.28) for i? = 0, we obtain the multiplicity of the diagram represe nted by M. This may be used to 
cross-check the multiplicities obtained before when solving the recursion relation (2.61). 



As an example, consider the following diagram of the four-point function with p — 3 vertices: 



(4.3) 



This diagram can be represented by altogether um — 3 different matrices, depending on the labeling of the top vertex 
with two external legs by 1, 2, or 3: 








1 


2 


3 







1 


2 


3 







1 


2 


3 








2 


1 


1 








1 


2 


1 








1 


1 


2 


1 


2 





1 


1 


1 


1 





1 


2 


1 


1 





2 


1 


2 


1 


1 





2 


2 


2 


1 





1 


2 


1 


2 





1 


3 


1 


1 


2 





3 


1 


2 


1 





3 


2 


1 


1 






(4.4) 



^From the zeroth row or column of these matrices or by inspecting the diagram ( |4.3| ), we read off that there exist two 
one-connections and one two-connection betwee n ex ternal legs and vertices. Thus Eq. ( [l.2| ) states that the number 
of identical vertex permutations of the diagram (13) is 4 (compare the corresponding entry in Table |l|). 



The matrix M contains of course all information on the topological properties of a diagram [^jo] . For this we define 
the submatrix M by removing the zeroth row and column from M . This allows us to recognize the connectedness of 
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a diagram: A diagram is disconnected if there is a vertex numbering for which M is a block matrix. Furthermore a 
vertex is a cutvertex, e.g. a vertex which hnks two otherwise disconnected parts of a diagram, if the matrix M has an 
almost block- form for an appropriate numbering of vertices in which the blocks overlap only on some diagonal element 
Mii, i.e. the matrix M takes block form if the ith row and column are removed. Similarly, the matrix M allows us 
to recognize a one-particle-reducible diagram, which falls into two pieces by cutting a certain line. Removing a line 
amounts to reducing the associated matrix elements in the submatrix M by one unit. If the resulting matrix M has 
block form, the diagram is one-particle-reducible. 



So far, the vertex numbering has been arbitrary, making the matrix representation of a diagram non-unique. To 
achieve uniqueness, we proceed as follows. We associate with each matrix a number whose digits are composed of the 
matrix elements M^j {0 < j < i < p), i.e. we form the number with the {p + l){p + 2)/2 — 1 elements 

MioMnAf2oM2iM22M3oM3iM32A/33 ...Mpp. (4.5) 



The smallest of these numbers is chosen to represent the diagram unique. For instance, the three matrices (4.4) carry 
the numbers 

201101120, 102101210, 101202110, (4.6) 

the smallest one being the last. Thus we uniquely represent the diagram (^^) by this number. 



B. Practical Generation 



We are now prepared for t he com puter generation of Feynman diagrams. First the vacuum diagrams are generated 
from the recursion relation (2.61). From these the diagrams of the connected two- and four-point functions are 
obtained by cutting or removing lines. A MATHEMATICA program perfor ms th is task. The resulting unique matrix 
representations of the diagrams up to the order p = 4 are listed in Table [vl- |VI]| . They are the same as those derived 



before by hand in Table |-[n. Higher-order results up to p = 6, containing all diagrams which are relevant for the 
five loop renormalization of the 0*-theory in d = 4 — e dimensions |^,|l^, are made available on internet |17|, where 
also the program can be found. 



1. Connected Vacuum Diagrams 



The computer solution of the recursion relation ( 2.61 ) necessitates to keep an exact record of the labeling of external 
legs of intermediate diagrams which arise from differentiating a vacuum diagram with respect to a line once or twice. 
To this end we have to extend our previous matrix representation of diagrams where the external legs are labeled 
as if they were connected to a simple additional vertex with number 0. For each matrix representing a diagram we 
define an associated vector which contains the labels of the external legs connected to each vertex. This vector has 
the length of the dimension of the matrix and will be prependcd to the matrix. By doing so, it is understood that 
the rows and columns of the matrix are labeled from to the number of vertices as ex pla ined in Subsection IV. A, so 
that we may omit these labels from now on. Consider, as an example, the diagram (4.3) of the four-point function 
with p = 3 vertices, where the spatial indices 1, 2, 3, 4 are assigned in a particular order: 




In our extended matrix notation, such a diagram can be represented in total by six matrices: 

/ {} 2 1 1\ / {} 1 2 1\ / {} 1 1 2\ 

{1,2} 2 1 1 {3} 10 12 {3} 10 2 1 

{3} 110 2' {1,2} 2 1 1 ' {4} 1201 

V {4} 112 0/ \ {4} 1210/ V {1,2} 2 1 1 0/ 

/ {} 2 1 1\ / {} 1 2 1\ / {} 112 

{1,2} 2 1 1 {4} 10 12 {4} 10 2 1 

{4} 110 2' {1,2} 2 1 1 ' {3} 1201 

V {3} 112 0/ V {3} 1210/ V {1,2} 2 1 1 
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I n th e calculation of the vacuum diagrams from the recursion relation (2.61), starting from the two-loop diagram 
( ^.62 ), we have to represent three different elementary operations in our extended matrix notation: 



1. Taking one or two derivatives of a vacuum diagram with respect to a line. For example, we apply this operation 
to the vacuum diagram #2 in Table | 



5Gi25G: 



34 



= 2 



<5Gi2 



-4 + 4 - 



(4.9) 



This has the matrix representation 



5Gi25G: 



34 











2 




3 


{1,3} 


2 









V{2,4} 


2 


2 






(4.10) 



The first and fourth matrix as well as the second and third matrix represent the same diagram in ( |4.9D , as can 
be seen bu permutating rows and columns of either matrix. 

2. Combining two or three diagrams to one. We perform this operation by creating a block matrix of internal lines 
from the submatrices representing the internal lines of the original diagrams. Then the zeroth row or column is 
added to represent the respective original external spatial arguments. Let us illustrate the combination of two 
diagrams by the example 



-o- 




{} 

{1,2} 



2 
2 1 




2 11 
2 10 
10 3 
10 3 



(4.11) 



and the combination of three diagrams by 



{} 
{1,2} 



2 
2 1 



{} 

{1,2,3,4} 



4 
4 



{} 
{3,4} 



2 
2 1 




2 4 2 

2 10 

4 

2 1 



(4.12) 



We observe that the ordering of the submatrices in the block matrix is arbitrary at this point, we just have to 
make sure to distribute the spatial labels of the external legs correctly. 

3. Connecting external legs with the same label and creating an internal line. This is achieved in our extended 
matrix notation by eliminating the spatial labels of external legs which appear twice, an d by performing an 
appropriate entry in the matrix for the additional line. Thus we obtain, for instance, from (4.11) 





111 

10 3 

13 



(4.13) 
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and similarily from ( 4.12| ) 



cxxo - 



{} 

U! 





12 

2 2 

2 1 



(4.14) 



As we reobtain at this stage connected vacuum diagrams where there are no more external legs to be labeled, 
we may omit the prepended vector. 

The selection of a unique matrix representation for the resulting vacuum diagrams obtained at each stage of the 
recursion relation proceeds as explained in detail in Subsection IV. A. By comparing we find out which of the vacuum 
diagr ams are topologically identical and sum up their individual multiplicities. Along these lines, the recursion rela- 
tion (2.61) is solved by a MATHEMATICA program up to the order p = 6. The results are shown in Table ^ and in 



Ref. [h7| . To each order p, the numbers of topologically different connected vacuum diagrams are 



p 


1 


2 


3 


4 


5 


6 


(0) 


1 


2 


4 


10 


28 


97 



(4.15) 



2. Two- and Four-Point Functions G12 and G1234 from Cutting Lines 



Having found all connected vacuum diagra ms, w e der ive from these the diagrams of the connected two- and four- 
point functions by using the relations (2.18) and ( ^.25 ). In the matrix representation, cutting a line is essentially 
identical to removing a line as explained above, except that we now interpret the labels which represent the external 
spatial labels as sitting on the end of lines. Since we are not going to distinguish between trivially "crossed" graphs 
which are related by exchanging external labels in our computer implementation, we need no longer carry around 
external spatial labels. Thus we omit the additional vector prepended to the matrix representing a diagram when 
generating vacuum diagrams. As an example, consider cutting a line in diagram ^3 in Table | 



SG- 



(4.16) 



which has the matrix representation 





(4.17) 



Here the plus signs and multiplication by 2 have a set-theoretical meaning and are not to be understood as matrix 
algebra operations. The last two matrices represent, incidentally, the same graph in ( 4.16 ) as can be seen by exchang- 
ing the last two rows and columns of either matrix. 

To create the connected four-point function, we also have to consider second derivatives of vacuum diagrams with 
respect to G~^. If an external line is cut, an additional external line will be created which is not connected to any 
vertex. It can be interpreted as a self-connection of the zeroth vertex which collects the external lines. This may be 
accomodated in the matrix notation by letting the matrix element Mqo count the number of lines not connected to 
any vertex. For example, taking the derivative of the first diagram in Eq. ( p^ ) gives 



Q 



SG- 



(4.18) 



with the matrix notation 



17 





(4.19) 



The first two matrices represent the same diagram as can be seen from Eq. ( 4.18 ). The last two matrices in Eq. ( 4.1£ ) 
correspond to disconnected diagrams: the first because of the absence of a connection between the two vertices, the 
second because of the disconnected line represented by the entry Mqo — 1. In the full expression for the two loop 
contribution 0^234 ^'^ four-point function in Eq. ( ^.25| ) all disconnected diagrams arising from cutting a line in 

(2) 

are canceled by diagrams resulting from the sum. Therefore we may omit the sum, take only the first term and 
discard all disconnected graphs it creates. This is particularly useful for treating low orders by hand. If we include 
the sum, we use the prescription of combining diagrams into one as described above in Subsection IV.B.l, except that 
we now omit the extra vector with the labels of spatial arguments. 



3. Two- and Four-Point Function G12 and G1234 from Removing Lines 

Instead of cutting lines of connected vacuum graphs once or twice, the pe rturb ati ve co efficie nts of G12 and G^234 
can also be obtained graphically by removing hues. Indeed, from ( ^.16D , ( ^.44D , ( p^ ) and ( ^.56D we get for the 
two-point function 



G12 — G12 



2 / G13G24 

'34 



so that we have for p > 



G 



12 



SG34 



2 / Gi3G24^t; — 

734 OU34 



(4.20) 



(4.21) 



at our disposal to compute the coefficients G^j'' from removing one line in the connected vacuum diagrams I^(P) in aU 
possible ways. The corresponding matrix operations are identical to the ones for cutting a line so that in this respect 
there is no difference between both procedures to obtain G12. 



Combining (4.2C) with (2.12), ( 2.23| ) and ( ^.56| ), we get for the connected four-point function 



^1234 — 4 / G15G26G37G48 J J 



'5678 

which is equivalent to 



G 



,ip) 

1234 



4 / Gi5G26G37G48 7-=; j-; — 

/5678 OLr5Q0Lr7s 



78 



p-1 

9=1 



4 / Gi5G27(G36G48 -I- G46G38)- 
"'5678 



SG56 SG 



78 



Gi5G27(G36G48 -|- G46G38) 



5678 



SGi 



56 



<5G78 



(4.22) 



(4.23) 



Again, the sum serves only to subtract disconnected diagrams which are created by the first term, so we may choose 
to discard both in the first term. 

Now the problem of generating diagrams is reduced to the generation of vacuum diagrams and subsequently taking 
functional derivatives with respect to G12. An advantage of this approach is that external lines do not appear at 
intermediate steps. So when one uses the cancellation of disconnected terms as a cross check, there are less operations 
to be performed than with cutting. At the end one just interprets external labels as sitting on external lines. Since 
all necessary operations on matrices have already been introduced, we omit examples here and just note that we can 
again omit external labels if we are not distinguishing between trivially "crossed" graphs. 

The generation of diagrams of the connected two- and four-point functions has been implemented in both possible 
ways. Cutting or removing one or two lines in the connected vacuum diagrams up to the order p — 6 lead to the 
following numbers n^p^ and n^^ of topologically different diagrams of G^2'' ^-i^^ ^1234- 



p 


1 


2 


3 


4 


5 


6 


(2) 
Up 


1 


3 


8 


30 


118 


548 


(4) 


1 


2 


8 


37 


181 


1010 



(4.24) 
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V. OUTLOOK 



In this work we have shown f that aU Feynman diagrams can be efficiently determined from a recursive graphical 
solution of a nonlinear functional differential equation order by order in the coupling strength. In separate publications, 
this method will be applied to the ordered phase of ^^-theory, where the energy functional contains a mixture of a 
cubic and a quartic interaction, and to quantum electrodynamics JTof . It is hoped that our method will eventually be 
combined with efficient numerical algorithms for actually evaluating these Feynman diagrams. 
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#2 
24 

(0,0,0,1;2) 



ceo 

(2,1,0,0;2) 



#4 
1728 

(0,3,0,0;6) 



#5 
3456 

(1,0,1,0;2) 



#6 
1728 

(3,0,0,0;6) 




(2,2,0,0;2) 



#8 
62208 

(0,4,0, 0;8) 




#9 
248832 

(0,2,0,0;8) 




#10 
55296 

(0,0,2,0;4) 



#11 

497664 

(1,2,0,0;2) 



#12 
165888 

(2,0,1,0;2) 



#13 
248832 

(2,1,0,0;4) 



#14 
165888 

(1,1,1,0;2) 



#15 
248832 

(3,1,0,0;2) 




#16 
62208 

(4,0,0,0;8) 




i?^'^'^ OXXX) 

(2,3,0,0;2) 



TABLE I. Connected vacuum diagrams and their multiplicities of the (/'''-theory up to five loops. Each diagram is char- 
acterized by the vector {S,D,T,F;N) whose components specify the number of self-, double, triple and fourfold connections, 
and of the identical vertex permutations, respectively. 
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#1.1 

12 Q 

(1,0, 0;2) 



#2.1 
192 

(0,0,1;2) 



#3.1 
288 

(1,1, 0;2) 



#3.2 
288 

(2,0,0;2) 



#4.1 
20736 

(0,2,0i2) 



#5.1 
6912 

(0,0, 1;4) 




#5.2 
20736 

(1,1, 0;2) 




#5.3 
13824 

(1,0,1;1) 



#6.1 
10368 

(2,0,0;4) 



#6.2 
10368 

(3,0, 0;2) 



#7.1 
10368 

(1,2,0;2) 



#7.2 
20736 

(2,1,0;1) 



#8.1 
995328 

(0,3,0;2) 



#9.1 

1990656 

(0,1, 0;4) 



#9.2 

1990656 

(0,2,0;2) 




#10.1 

221184 

(0,0, 2;2) 



#10.2 

663552 

(0,1, 1;2) 



#11.1 

995328 

(0,2,0;4) 



#11.2 

1990656 Q_ 

(1,2,0;1) 



#11.3 
995328 

(1,2,0;2) 




#11.4 
3981312 

(1,1,0;1) 



#12.1 

995328 

(2,1,0;2) 



#12.2 
331776 

(2,0,1;2) 



#12.3 
663552 

(2,0, 1;1) 



#12.4 
663552 

(1,0,1;2) 




#13.1 
995328 

(2,0,0;4) 



#13.2 
995328 

(1,1, 0;4) 



#13.3 
1990656 

(2,1,0;1) 



#14.1 
995328 

(1,2,0;2) 




#14.2 
663552 

(i,i,i;i) 



■e- 



#14.3 
663552 

(1,0,1;2) 



a 



#14.4 
331776 

(0,1, 1;4) 



#15.1 
995328 

(3,1,0;1) 




#15.2 
497664 

(3,1,0;2) 




#15.3 
497664 

(2,1,0;4) 



#15.4 
995328 

(2,1,0;2) 



#15.5 
995328 

(3,0,0;2) 



#16.1 
497664 

(3,0,0;4) 




#16.2 
497664 

(4,0,0i2) 



#17.1 
497664 

(1,3,0;2) 



#17.2 
995328 

(2,2,0;1) 




#17.3 
497664 

(2,2,0i2) 




TABLE II. Connected diagrams of the two-point function and their multiplicities of the (/>'*-theory up to four loops. Each 
diagram is characterized by the vector {S,D,T;N) whose components specify the number of self-, double, triple connections, 
and of the identical vertex permutations, respectively. 



21 



■yC,(p) 



#1.1.1 

24 

(0,0,0;24) 



#2.1.1,#3.1.1 
1152,576 

1728 
(0,1, 0;8) 



XX 



#3. 1.2, #3. 2.1 
1162,1152 

2304 
{1,0, 0;6) 



#4.1.1, #7.1.1 

41472,20736 
62208 
(0,2,0i8) 



#4.1. 2, #5. 1.1, #5. 2.1 
165888,41472,41472 
248832 
(0,1, 0;4) 



#5. 1.2, #5. 3. 2 
27648,27648 
55296 
(0,0, 1;6) 



#5. 2. 2, #6.1.1 
82944,41472 
124416 
(1,0, 0;8) 



#5. 2. 3, #5. 3.1, #7. 1.2, #7. 2.1 
82944,82944,41472,41472 
248832 
(1,1, 0;2) 



#6. 1.2, #6. 2. 2, #7. 2. 2 
20736,20736,82944 
124416 
(2,0, 0;4) 



QiQ 



#6. 1.3, #6.2.1 
41472,41472 

82944 
(2,0,0;6) 



#7.1. 3, #7.2. 3 
41472,41472 

82944 
(1,1, 0;6) 



#8. 1.1, #17.1.1 
1990656,995328 
2985984 
(0,3,0;8) 



#8. 1.2, #9. 2.1, #10. 2.1 

3981312,3981312,3981312 
11943936 
(0,2,0;4) 



#8. 1.3, #11. 1.2, #11. 3.1 

7962624,1990656,1990656 
11943936 
(0,2,0;4) 



#9. 1.1, #13. 2.1 
3981312,1990656 
5971968 
(0,1,0;16) 




#9.1.2 

7962624 
(0,0,0;24) 




#9. 1.3, #9. 2. 3, #11. 1.1, #11. 4.1 
15925248,15925248,7962624,7962624 
47775744 
(0,1, 0;2) 




#9. 2. 2, #14. 1.1, #14.4.3 
7962624,1990656,1990656 

11943936 
(0,2,0;4) 



#10. 1.1, #10. 2. 3, #14.2.1, #14.4.2 
2654208,2654208,1327104,1327104 

7962624 
(0,1, 1;2) 



#10. 2. 2, #12. 4.1 
2654208,1327104 

3981312 
(0,0, 1;8) 



#11. 1.3, #11. 2.1 

3981312,3981312 
7962624 
(0,2,0;6) 



#11. 2. 2, #11. 4. 3, #14. 1.2, #14.3. 3 
7962624,7962624,3981312,3981312 

23887872 

(1,1, 0;2) 




#11. 2. 3, #11.4.2, #13. 2. 2, #13. 3.1 
7962624,7962624,3981312,3981312 

23887872 

(1,1, 0;2) 
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#11. 2. 4,#11. 3. 2, #17. 1.2, #17.2.1 
3981312,3981312,1990656,1990656 

11943936 

(1,2,0;2) 



#11. 4. 6, #13. 1.1, #13. 2. 3 
15925248,3981312,3981312 
23887872 
(1,0,0;4) 



#12.1.4,#12.3.3,#15.1.1,#16.3.2 
3981312,3981312,1990656,1990656 

11943936 
(2,1,0;2) 



#12. 3.1, #12. 4.4 

1327104.1327104 
2654208 
(1,0, 1;6) 





#13. 2. 4, #13. 3. 5 
3981312,3981312 
7962624 
(1,1, 0;6) 




#11. 3. 3, #11.4.4, #12. 1.1, #12.4.6 
7962624,7962624,3981312,3981312 

23887872 

(1,1, 0;2) 



#12. 1.2, #12. 2. 2, #13. 3. 3, #17. 2. 2 
1990656,1990656,3981312,3981312 

11943936 

(2,1,0;2) 



#12. 2.1, #12.4.2 
1327104,1327104 

2654208 
(1,0, 1;6) 



#13. 1.2, #13. 3. 4, #15. 4.2, #15. 6.1 
7962624,7962624.3981312,3981312 

23887872 

(2,0, 0;2) 



#13. 3. 2, #15. 2.1, #15. 3. 3 
3981312,995328,995328 
5971968 
(2,1,0;4) 




#11.4.6, #15. 3.1, #15.4.1 
7962624,1990656,1990656 

11943936 

(1,1, 0;4) 



#12. 1.3, #16. 1.2 
3981312,1990656 

5971968 

(2,0,0;8) 



#12. 3. 2, #12. 4. 3, #14. 2. 2, #14.3. 2 
1327104,1327104,2654208,2654208 

7962624 

(1,0, 1;2) 



#13. 1.3, #16. 1.1 
1990656,995328 
2985984 
(2,0,0;16) 



#14.1.3,#14.2.3,#17.1.:5,#17.3.1 
3981312,3981312,1990656,1990656 

11943936 

(1,2,0;2) 





#14.1. 4,#15. 4.4 
3981312,1990656 

5971968 

(1,1, 0;8) 




#14. 3.1, #14.4.1 
1327104,1327104 

2654208 
(0,0,1;12) 



#15. 1.2, #15. 5. 3, #16. 1.3, #16. 2. 2 
3981312,3981312,1990656,1990656 

11943936 
(3,0,0;2) 



#15. 1.3, #15. 4. 3, #17.2. 3, #17.3. 2 
1990656,1990656,3981312,3981312 

11943936 

(2,1,0;2) 



#15.2.3,#18.4.6 
1990656,1990656 

3981312 

(2,1,0;6) 



#15. 1.4, #15. 4. 5 
1990656,1990656 

3981312 

(2,1,0;6) 



#15. 2. 2, #15. 5. 2 

1990656,1990656 
3981312 
(3,0,0;6) 



#15. 3. 4, #15. 5. 4 
1990656,1990656 
3981312 
(2,0,0; 12) 



#16. 1.4, #16. 2.1 
1990656,1990656 

3981312 

(3,0,0;6) 



#17. 1.4, #17. 2. 4 
1990656,1990656 

3981312 

(1,2,0;6) 




+0- 





TABLE III. Connected diagrams of the four-point function and their multiplicities of the c'i^-theory up to three loops. Each 
diagram is characterized by the vector {S, D, T; N) whose components specify the number of self-, double, triple connections, 
and of the identical vertex permutations, respectively. 



23 



#1.1 

2 



1 oo 



#2.1 
48 

(0,0;8) 



#2.2 
24 

(0,2;4) 



OD 



#3.1 
96 

(0,0;4) 



#3.2 
96 

(1,0;2) 



(XK) 



#3.3 

24 
(2,i;2) 



CKX) 



#4.1 
3840 

(0,0;12) 

#5.1 

11520 

(0,1;2) 



#4.2 
11520 

(0,0i4) 



#5.2 

23040 

(0,0;2) 



#4.3 
3840 

(0,0; 12) 



#4.4 

(0,3;6) 



#5.3 

11520 

(1,0;2) 



#5.4 

5760 

(l,li2) 



#4.5 
5760 

(0,1;4) 



CHID 



#6.1 

7680 

(0,0;6) 



#6.2 
11520 

(1,0;2) 



gK3 



#6.3 
5760 

(2,0;2) 



#6.4 
960 

(3,0;6) 




#7.1 
11520 

(0,0;4) 



#7.2 
5760 

(0,0;8) 



(M) 



#7.3 
11520 

(1,0;2) 



QK3 



#7.4 
5760 

(i,i;2) 



(IKX3 



#7.5 
1440 

(2,2;2) 



#7.6 
2880 

(2,0;4) 



TABLE IV. Connected vacuum diagrams and their multiplicities of the 0^ ^-theory up to four loops. Each diagram is 
characterized by the vector (S, D; N) whose components specify the number of self- and double connections as well as the 
identical vertex permutations, respectively. 



24 



W^^^: 1 diagram 


i 


1 








j 


1 








# 




{S,D,T,F;N) 


M 




1 


2 


(2,1,0,0;1) 


3 


8 



W'''^^: 2 diagrams 


i 


1 


22 








j 


1 


12 








# 


Mij 


{S,D,T,F;N) 


M 




2 





40 


(0,0,0,1;2) 


24 


48 


3 


1 


21 


(2,1,0,0;2) 


72 


16 



W^^': 4 diagrams 


i 


1 


22 


333 








j 


1 


12 


123 








# 




(S,D,T,F;N) 


M 




5 





11 


310 


(1,0,1,0;2) 


3456 


24 


4 





20 


220 


(0,3,0,0;6) 


1728 


48 


7 





21 


201 


(2,2,0,0;2) 


2592 


32 


6 


1 


11 


111 


(3,0,0,0;6) 


1728 


48 



W''*': 10 diagrams 



i 


1 


22 


333 


4444 








j 


1 


12 


123 


1234 








# 




{S,D,T,F;N) 


M 




10 





00 


130 


3100 


(0,0,2,0;4) 


55296 


144 


8 





00 


220 


2200 


(0,4,0,0;8) 


62208 


128 


12 





01 


111 


3100 


(2,0,1,0;2) 


165888 


48 


14 





01 


120 


3010 


(1,1,1,0;2) 


165888 


48 


17 





01 


201 


2200 


(2,3,0,0;2) 


124416 


64 


11 





01 


210 


2110 


(1,2,0,0;2) 


497664 


16 


9 





10 


120 


2110 


(0,2,0,0;8) 


248832 


32 


13 





11 


101 


2110 


(2,1,0,0;4) 


248832 


32 


15 





11 


111 


2001 


(3,1,0,0;2) 


248832 


32 


16 


1 


01 


111 


1101 


(4,0,0,0;8) 


62208 


128 



TABLE V. Unique matrix representation of all connected vacuum diagrams of <^*-theory up to the order p — 4. The number 
in the first column corresponds to their graphical representation in Table ^. The matrix elements Mij represent the numbers of 
lines connecting two vertices i and j, with omitting Mio — for simplicity. The running numbers of the vertices are listed on 
top of each column in the first two rows. The further columns contain the vector {S, D, T, F; N) characterizing the topology of 
the diagram, the multiplicity M and the weight W — Af/[(4!)''p!]. 



25 



Gji,"*: 1 diagram 


i 


11 








j 


01 








# 




{S,D,T;N) 


M 




1.1 


21 


(1,0,0;2) 


12 


4 







'^12 


: 3 diagrams 




i 


11 


222 








j 


01 


012 








# 


Mij 


(S,D.T:N) 


M 




3.1 


01 


220 


(1,1,0;2) 


288 


8 


2.1 


10 


130 


(0,0,1;2) 


192 


12 


3.2 


11 


111 


(2,0,0;2) 


288 


8 







•^12 • 


8 diagrams 




i 


11 


222 


3333 








j 


01 


012 


0123 








# 




(S,D,T;N) 


M 




7.1 


00 


021 


2200 


(1,2,0;2) 


10368 


16 


5.1 


00 


030 


2110 


(0,0,1;4) 


6912 


24 


5.3 


00 


111 


1300 


(1,0,1;1) 


13824 


12 


4.1 


00 


120 


1210 


(0,2,0;2) 


20736 


8 


6.1 


01 


Oil 


2110 


(2,0,0;4) 


10368 


16 


7.2 


01 


101 


1210 


(2,1,0;1) 


20736 


8 


5.2 


01 


110 


1120 


(1,1,0;2) 


20736 


8 


6.2 


01 


111 


1101 


(3,0,0;2) 


10368 


16 



Gjj' '• 30 diagrams 



i 


11 


222 


3333 


44444 










j 


01 


012 


0123 


01234 










# 




(S,D,T 


JV) 


M 




17.1 


00 


001 


0220 


22000 


(1,3,0 


2) 


497664 


32 


12.4 


00 


001 


0310 


21100 


(1,0,1 


2) 


663552 


24 


14.2 


00 


001 


1120 


13000 


(1,1,1 


1) 


663552 


24 


11.3 


00 


001 


1210 


12100 


(1,2,0 


2) 


995328 


16 


14.4 


00 


010 


0130 


22000 


(0,1,1 


4) 


331776 


48 


11.1 


00 


010 


0220 


21100 


(0,2,0 


4) 


995328 


16 


10.1 


00 


010 


1030 


13000 


(0,0,2 


2) 


221184 


72 


9.2 


00 


010 


1120 


12100 


(0,2,0 


2) 


1990656 


8 


15.3 


00 


on 


0111 


22000 


(2,1,0 


4) 


497664 


32 


15.4 


00 


on 


0201 


21100 


(2,1,0 


2) 


995328 


16 


13.2 


00 


on 


0210 


21010 


(1,1,0 


4) 


995328 


16 


12.3 


00 


on 


1011 


13000 


(2,0,1 


1) 


663552 


24 


13.3 


00 


on 


1101 


12100 


(2,1,0 


1) 


1990656 


8 


11.4 


00 


on 


1110 


12010 


(1,1,0 


1) 


3981312 


4 


11.2 


00 


020 


1011 


12100 


(1,2,0 


1) 


1990656 


8 


8.1 


00 


020 


1020 


12010 


(0,3,0 


2) 


995328 


16 


9.1 


00 


020 


1110 


11110 


(0,1,0 


4) 


1990656 


8 


17.2 


00 


021 


1001 


12010 


(2,2,0 


1) 


995328 


16 


14.1 


00 


021 


1100 


11020 


(1,2,0 


2) 


995328 


16 


15.2 


00 


021 


1101 


11001 


(3,1,0 


2) 


497664 


32 


14.3 


00 


030 


1001 


11110 


(1,0,1 


2) 


663552 


24 


10.2 


00 


030 


1010 


11020 


(0,1,1 


2) 


663552 


24 


12.2 


00 


030 


1011 


11001 


(2,0,1 


2) 


331776 


48 


16.1 


01 


001 


0111 


21100 


(3,0,0 


4) 


497664 


32 


15.1 


01 


001 


1011 


12100 


(3,1,0 


1) 


995328 


16 


17.3 


01 


001 


1020 


12010 


(2,2,0 


2) 


497664 


32 


13.1 


01 


001 


1110 


11110 


(2,0,0 


4) 


995328 


16 


15.5 


01 


on 


1001 


11110 


(3,0,0 


2) 


995328 


16 


12.1 


01 


on 


1010 


11020 


(2,1,0 


2) 


995328 


16 


16.2 


01 


on 


1011 


11001 


(4,0,0 


2) 


497664 


32 



TABLE VI. Unique matrix representation of all connected two-point function of (/)''-theory up to the order p = 4. The 
numbers in the first column correspond to their graphical representation in Table ^ The matrix elements Alij represent the 
numbers of lines connecting two vertices i and j, with omitting Mio = for simplicity. The running numbers of the vertices 
are listed on top of each column in the first two rows. The further columns contain the vector {S, D, T; N) characterizing the 
topology of the diagram, the multiplicity M and the weight W = M /[{4\Yp\]. 



26 



G^234 • 1 diagram 


i 


11 








j 


01 








# 




(S,D,T;N) 


M 




1.1.1 


40 


(0,0,0;24) 


24 


24 



^1234- 2 diagrams 


i 


11 


222 








j 


01 


012 








# 




(S,D.T:N) 


M 




3.1.2, 3.2.1 


11 


310 


(1,0,0;6) 


2304 


12 


2.1.1, 3.1.1 


20 


220 


(0,1,0;8) 


1728 


16 



"-'1234 • 


8 d 


lagrams 






1 


11 


222 


3333 








j 


01 


012 


0123 








-U. 

W 


M.J 


{S,D ,T;N) 


M 


w-^ 


CC 1 O CC Q O 

o.i.z, 0.O.2 


nn 
uu 




1 nn 

OiUU 


(0,U,l,bj 


55296 


36 


/111 VII 

4.1.1, /.l.l 


00 


220 


2200 


(,0,2,0,»j 


62208 


32 


D.l.o, D.Z.I 


01 


111 


3100 


(,2,U,0,0j 


82944 


24 


7.1.3, 7.2.3 


01 


120 


3010 




82944 


24 


5.2.2, 6.1.1 


01 


210 


2110 


(1,0,0;8) 


124416 


16 


5.2.3, 5.3.1, 7.1.2, 7.2.1 


10 


111 


2200 


(1,1,0;2) 


248832 


8 


4.1.2, 5.1.1, 5.2.1 


10 


120 


2110 


(0,1,0;4) 


248832 


8 


6.1.2, 6.2.2, 7.2.2 


11 


101 


2110 


(2,0,0;4) 


124416 


16 



G 



37 diagrams 





i 


11 


222 


3333 


44444 










3 


01 


012 


0123 


01234 








# 


Mij 


(S.D,T;N) 


M 






13.2.4, 13.3.5 


00 


oil 


1210 


31000 


(1,1,0;6) 


7962624 


24 




12.3.1, 12.4.4 


00 


oil 


1300 


30100 


(1,0,1;6) 


2654208 


72 




11.4.5, 15.3.1, 15.4.1 


00 


oil 


2110 


22000 


(1,1,0;4) 


11943936 


16 




11.1.3, 11.2.1 


00 


020 


1120 


31000 


(0,2,0;6) 


7962624 


24 




8.1.1, 17.1.1 


00 


020 


2020 


22000 


(0,3,0;8) 


2985984 


64 




9.1.1, 13.2.1 


00 


020 


2110 


21100 


(0,1,0;16) 


5971968 


32 




15.2.3, 15.4.6 


00 


021 


1101 


31000 


(2,1,0;6) 


3981312 


48 




17.1.4, 17.2.4 


00 


021 


1200 


30010 


(1,2,0;6) 


3981312 


48 




14.1.4, 15.4.4 


00 


021 


2100 


21010 


(1,1,0;8) 


5971968 


32 




12.2.1, 12.4.2 


00 


030 


1011 


31000 


(1,0,1;6) 


2654208 


72 




14.3.1, 14.4.1 


00 


030 


1110 


30010 


(0,0,1;12) 


2654208 


72 




10.2.2, 12.4.1 


00 


030 


2010 


21010 


(0,0,1;8) 


3981312 


48 


11 


2.4, 11.3.2, 17.1.2, 17.2.1 


00 


101 


1210 


22000 


(1,2,0;2) 


11943936 


16 


12 


3.2, 12.4.3, 14.2.2, 14.3.2 


00 


101 


1300 


21100 


(1,0,1;2) 


7962624 


24 




9.2.2, 14.1.1, 14.4.3 


00 


110 


1120 


22000 


(0,2,0;4) 


11943936 


16 


9 


1.3, 9.2.3, 11.1.1, 11.4.1 


00 


110 


1210 


21100 


(0,1,0;2) 


47775744 


4 


10 


1.1, 10.2.3, 14.2.1, 14.4.2 


00 


110 


1300 


20200 


(0,1,1;2) 


7962624 


24 




13.3.2, 15.2.1, 15.3.3 


00 


111 


1101 


22000 


(2,1,0;4) 


5971968 


32 


11 


2.3, 11.4.2, 13.2.2, 13.3.1 


00 


111 


1200 


21010 


(1,1,0;2) 


23887872 


8 




8.1.3, 11.1.2, 11.3.1 


00 


120 


1200 


20110 


(0,2,0;4) 


11943936 


16 




15.1.4, 15.4.5 


01 


001 


1120 


31000 


(2,1,0;6) 


3981312 


48 




13.1.3, 16.1.1 


01 


001 


2110 


21100 


(2,0,0;16) 


2985984 


64 




16.1.4, 16.2.1 


01 


Oil 


1011 


31000 


(3,0,0;6) 


3981312 


48 




15.3.4, 15.5.4 


01 


Oil 


1110 


30010 


(2,0,0;12) 


3981312 


48 




12.1.3, 16.1.2 


01 


Oil 


2010 


21010 


(2,0,0;8) 


5971968 


32 


11 


3.3, 11.4.4, 12.1.1, 12.4.5 


01 


100 


1120 


21100 


(1,1,0;2) 


23887872 


8 


14 


1.3, 14.2.3, 17.1.3, 17.3.1 


01 


100 


1210 


20200 


(1,2,0;2) 


11943936 


16 


15 


1.2, 15.5.3, 16.1.3, 16.2.2 


01 


101 


1101 


21100 


(3,0,0;2) 


11943936 


16 


13 


1.2, 13.3.4, 15.4.2, 15.5.1 


01 


101 


1110 


21010 


(2,0,0;2) 


23887872 


8 


15 


1.3, 15.4.3, 17.2.3, 17.3.2 


01 


101 


1200 


20110 


(2,1,0;2) 


11943936 


16 


12 


1.4, 12.3.3, 15.1.1, 15.3.2 


01 


110 


1101 


20200 


(2,1,0;2) 


11943936 


16 




11.4.6, 13.1.1, 13.2.3 


01 


110 


1110 


20110 


(1,0,0;4) 


23887872 


8 




8.1.2, 9.2.1, 10.2.1 


10 


100 


1120 


12100 


(0,2,0;4) 


11943936 


16 


12 


1.2, 12.2.2, 13.3.3, 17.2.2 


10 


101 


1101 


12100 


(2,1,0;2) 


11943936 


16 


11 


2.2, 11.4.3, 14.1.2, 14.3.3 


10 


101 


1110 


12010 


(1,1,0;2) 


23887872 


8 




9.1.2 


10 


110 


1110 


11110 


(0,0,0;24) 


7962624 


24 




15.2.2, 15.5.2 


10 


111 


1101 


11001 


(3,0,0;6) 


3981312 


48 



TABLE VII. Unique matrix representation of all connected two-point function of 0^-theory up to the order p = 4. The 
numbers in the first column correspond to their graphical representation in Table |l|. The matrix elements Mij represent the 
numbers of lines connecting two vertices i and j, with omitting Mio = for simplicity. The running numbers of the vertices 
are listed on top of each column in the first two rows. The further columns contain the vector {S, D, T; N) characterizing the 
topology of the diagram, the muItipHcity M and the weight W = M /[{4\)^p\]. 
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